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We consider the linear stability of dissipative MHD Taylor-Couette flow with imposed toroidal 
magnetic fields. The inner and outer cylinders can be either insulating or conducting; the inner 
one rotates, the outer one is stationary. The magnetic Prandtl number can be as small as 10~^, 
approaching realistic liquid-metal values. The magnetic field destabilizes the fiow, except for radial 
profiles of B^{R) close to the current-free solution. The profile with Bin ~ Bout (the most uniform 
field) is considered in detail. For weak fields the TC-fiow is stabilized, until for moderately strong 
fields the m = 1 azimuthal mode dramatically destabilizes the flow again. There is thus a maximum 
value for the critical Reynolds number. For sufficiently strong fields (as measured by the Hartmann 
number) the toroidal field is always unstable, even for Re = 0. 

The electric currents needed to generate the required toroidal flelds in laboratory experiments are 
a few kA if liquid sodium is used, somewhat more if gallium is used. Weaker currents are needed for 
wider gaps, so a wide-gap apparatus could succeed even with gallium. The critical Reynolds numbers 
are only somewhat larger than the nonmagnetic values, so such an experiment would require only 
modest rotation rates. 

PACS numbers: 47.20.Ft, 47.65. -fa 



I. MOTIVATION 

Taylor-Couette flows of electrically conducting fluid 
between rotating concentric cylinders are a classical prob- 
lem of hydrodynamic and hydromagnetic stability the- 
ory. It is becoming increasingly clear that the stability 
of differential rotation combined with magnetic fields is 
also one of the key problems in MHD astrophysics. For 
uniform axial magnetic fields this phenomenon is called 
the magnetorotational instability (MRI). The Keplerian 
rotation of accretion disks with weak vertical magnetic 
fields becomes unstable, so that angular momentum is 
transported outward (star formation) and gravitational 
energy is efficiently transformed into heat and radia- 
tion (quasars) . Galactic rotation profiles (approximately 
51 cx may also be MRI-unstable, resulting in inter- 

stellar MHD turbulence [l|,[l|. 

Due to the small size of laboratory experiments, and 
the extremely small magnetic Prandtl numbers of liquid 
metals, it is very difficult to achieve the MRI in the lab 
with a purely axial magnetic field. However, by adding a 
toroidal field i?^ cx (which is current-free within the 
fluid), it was possible to obtain the MRI experimentally 
0) 01 ) with the most unstable mode being axisymmetric 
and oscillatory, just as predicted theoretically Q. 
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If the field strength is sufficiently great, such a current- 
free toroidal field yields a magnetorotational instability 
even without an axial field; the instability is simply non- 
axisymmctric rather than axisymmetric. The simulta- 
neous occurrence of a stable differential rotation and a 
stable toroidal field may therefore nevertheless be un- 
stable, provided that the magnetic Reynolds number ex- 
ceeds O(IOO) [i]. 

In contrast to these magnetorotational instabilities, in 
which the magnetic field acts as a catalyst, but not as a 
source of energy, toroidal fields that are not current-free 
may become unstable more directly, by so-called current 
or Tayler instabilities 0, Q • Because the source of en- 
ergy is now the current rather than the differential ro- 
tation, these (non-axisymmetric) instabilities can exist 
even without any differential rotation, provided only the 
current is large enough. The magnitude depends quite 
strongly on the precise radial profile of the associated 
magnetic field B^{R), but not on the magnetic Prandtl 
number. The topic of this paper is how Tayler instabil- 
ities interact with differential rotation, and whether it 
might be possible to realize some of the resulting modes 
in laboratory experiments. The combination of Tayler in- 
stabilities and differential rotation may also be relevant 
to a broad range of astrophysical problems, including the 
stability of the solar tachocline, the existence of active 
solar longitudes Q, the flip-flop phenomenon of stellar 
activity [lO|, A-star magnetism [If], and even the possi- 
bility of 'a differential rotation driven dynamo in a stably 
stratified star' jl^ . 
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II. EQUATIONS 

According to the Rayleigh criterion an ideal flow is 
stable against axisymmetric perturbations whenever the 
specific angular momentum increases outwards 



(1) 



where (R, 0, z) are cylindrical coordinates, and is the 
angular velocity. The necessary and sufficient condition 
for the axisymmetric stability of an ideal Taylor-Couette 
flow with an imposed azimuthal magnetic fleld Ba, is 



R^dR^ ' fioP dR\R 



>0, 



(2) 



where is the permeability and p the density p^ . [3] ■ 
In particular, all ideal flows can thus be destabilized, by 
azimuthal magnetic fields with the right profiles and am- 
plitudes. Any fields increasing outward more slowly than 
ex R, including in particular the outwardly decreasing 
current-free field B^ cx have a stabilizing infiuence 
though [ii . 

Tayler Q found the necessary and sufficient condition 



dR 



{RBl) > 0. 



(3) 



for the non-axisymmetric stability of an ideal fluid at 
rest. Outwardly increasing fields are therefore unstable 
now (but -B^ oc is still stable). If this condition 
([3]) is violated, the most unstable mode has azimuthal 
wavenumber m = 1. In this paper we wish to consider 
how these Tayler instabilities are modified if the fluid is 
not at rest, but is instead differentially rotating. 

We will find that, depending on the magnitudes of 
the imposed differential rotation and magnetic fields, and 
also on the magnetic Prandtl number, a magnetic field 
may either stabilize or destabilize the differential rota- 
tion, and the most unstable mode may be either the ax- 
isymmetric Taylor vortex flow, or the non-axisymmetric 
Tayler instability. We focus on the limit of small mag- 
netic Prandtl numbers appropriate for liquid metals, and 
calculate the rotation rates and electric currents that 
would be required to obtain some of these instabilities 
in liquid metal laboratory experiments. 

The governing equations are 

^ + (UV)U = --VP + vAU + —curl B x B, (4) 
ot p Ho 



5B 

'dt 



curl(U X B) -I- ryAB, 



and 



div U = div B = 0, 



(5) 



(6) 



where U is the velocity, B the magnetic field, P the pres- 
sure, V the kinematic viscosity, and 77 the magnetic dif- 
fusivity. 

The basic state \s Ub = Uz = Bp = B^ = and 



U^ = Rn = anR 



R 



B^^aBR+^, (7) 



where ao, 60, as and are constants defined by 

an — "in—; On = "in-Kjjj^^ TJ, 

1 - flBV 



as 



where 
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R'm ^ ^out ^ -^out 



Ro 



, (8) 



(9) 



Rin and i?out are the radii of the inner and outer cylin- 
ders, r^in and r^out are their rotation rates (we will in fact 
fix r^out = for all results presented here), and Bm and 
Bout are the azimuthal magnetic fields at the inner and 
outer cylinders. The possible magnetic field solutions are 
plotted in Fig. [1] Note though that - unlike fi, where 
r^in and Qout are the physically relevant quantities - for 
B^ the fundamental quantities are not so much Bin and 
-Bout, but rather as and bs themselves. In particular, a 
field of the form bs/R is generated by running an axial 
current only through the inner region R < i?i„, whereas 
a field of the form asR is generated by running an axial 
current through the entire region R < Rout , including the 
fluid. One of the aspects we will be interested in later on 
is how large these currents must be, and whether they 
could be generated in a laboratory experiment. 

We are interested in the linear stability of the basic 
state ([7]). The perturbed quantities of the system are 
given by 



ur, RVl + UA, Uz, bR, Bs + b^, 62 



(10) 



Applying the usual normal mode analysis, we look for 
solutions of the linearized equations of the form 



F = F(i?)exp(i(A:z + mcj) + ujt)) . 



(11) 



The dimensionless numbers of the problem are the mag- 
netic Prandtl number Pm, the Hartmann number Ha, 
and the Reynolds number Re, given by 



Pm = — , 
77 



Ha = 



BinRo 



Re = 



(12) 



where Rq = (RiniRout — Rin)y^^ is the unit of length. 

Using (|lip . linearizing the equations Q and ([5]), and 
representing the result as a system of first order equa- 
tions, we have 

duR ur .m 
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where X2, and are defined as 
X2 



du^ _duz _db^ b^ 

-dR^li' ^'-dR^R-^^^^ 



Length has been scaled by Rq, time by ilj^"'^, the basic 
state angular velocity by fiin, the perturbation velocity 
by rj/Ro, and the magnetic fields, both basic state and 
perturbation, by Bin- 
Am appropriate set of ten boundary conditions is 
needed to solve the system p^ . For the velocity the 
boundary conditions are always no-slip. 



Ur = Ua, ^ Uz ^ 0. 



(15) 



For the magnetic field the boundary conditions depend 
on whether the walls are insulators or conductors. For 
conducting walls the radial component of the field and 
the tangential components of the current must vanish, 
yielding 



db^/dR + ^bR^Q. 



(16) 



These boundary conditions are applied at both i?in and 

For insulating walls the boundary conditions are some- 
what more complicated; matching to interior and exterior 
potential fields then yields 



TO 

kR 



(17) 
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/„(fci?)+/™+i(fci?)) =0, (18) 



(19) 
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JR ■ 



Kra{kR) 



(^if™(fci?) - Krn+l{kR)') = (20) 



at i? = i?outj where /„ and Kn are the modified Bessel 
functions ui 
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FIG. 1: The basic state profiles of B^{R), for 77 = 0.5. Re- 
membering that -B^(l) has been normalized to 1, we find that 
fiB is given simply by B^{2), and can therefore be read off the 
right-hand axis. The stability domain for m = 1 (see (I2ip ') 
is cross-hatched; the cross-hatched and hatched domains to- 
gether are the stability domain for m = (see (HU). The 
current-free solution — 1/R is given by the gray line 
{fiB = 0.5). The electric currents inside and outside the inner 
cylinder are parallel above (Ib ~ 0.5 and anti-parallel below 
ftB = 0.5, that is, the signs of ub and 6_b are the same above 
0.5, and opposite below. 

Given the basic state ([7]), Tayler's stability condition 
to nonaxisymmetric perturbations becomes 



< As < 



477(1 



3 — 27)2 — fj 



fi4 



Ml- 



(21) 



Note that /ti 1 (but is always less than 1) if — > 1. 
For 77 = 0.5 we have — 0.62. Similarly, the stability 
condition to axisymmetric perturbations becomes 



< As < 4 = Mo- 



(22) 



For fj — 0.5 we have /to — 2. For < 77 < 1 we always 
have fii < fio, so that the stability interval pT|) for to = 1 
is much smaller than the stability interval (|22|) for to = 0, 
as shown in Fig. [TJ The current-free solution fiB — 0.5 
is of course always stable. 
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III. BASIC RESULTS 

Figure [2] shows how the stability curves depend on Ha 
and Re, for Pm = 10, 1, 0.1 and 0.01, with Ao = 
(stationary outer cylinder) and /is = 1 {B^ as uniform 
as possible). For Ha = we find that m = goes un- 
stable before m = 1, at the critical Reynolds number 
Recrit = 68; this is just the familiar value for the on- 
set of nonmagnetic Taylor vortices (at this particular ra- 
dius ratio). Being entirely nonmagnetic, this value obvi- 
ously does not depend on Pm. At the other limiting case, 
Re — 0, we find that only m — 1 goes unstable, at the 
critical Hartmann number Hacrit = 150. These Tayler in- 
stabilities also turn out to be independent of Pm, despite 
being driven by the magnetic field. 

We are interested in how these two limiting cases Ha = 
and Re = are connected, and how the two types 
of instabilities interact when neither parameter is zero. 
The two instabilities certainly are connected; the m = 1 
modes in the two limiting cases are smoothly joined to 
one another for all Prandtl numbers. The nature of the 
interaction is quite different though, depending on Pm. 

Turning to Pm = 1 first, we see that there is relatively 
little interaction between rotational and magnetic effects; 
instability simply sets in as long as either Ha > Hacrit or 
Re > Rccrit. For Pm = 10 the situation is very dif- 
ferent. There we find a broad range of parameters, for 
example Ha = 100 and Re — 50, that would be stable 
if rotational or magnetic effects were acting alone, but 
which are now unstable, due to the interaction between 
the two (see also 17]). Finally, for Pm = 0.1 we have 
the opposite situation, namely a range of parameters, for 
example Ha — 100 and Re — 100, that would be unstable 
if rotational effects were acting alone, but which arc now 
stable. 

Small Pm are generally stabilizing. The opposite is 
true for Pm > 1. As shown in Fig. [2] (top) instability 
then also sets in for Hartmann numbers less than 150. In 
other words, for Hartmann numbers exceeding around 50, 
the critical Reynolds number for the onset of instability 
is much smaller than 68. 



IV. LIQUID METALS 

Having explored the general behavior for a range of 
magnetic Prandtl numbers, we now focus attention on 
the limit of very small Pm, such as would apply for exper- 
iments involving liquid metals. We will here consider con- 
ducting and insulating boundary conditions separately. 



A. Conducting cylinder walls 

Figure [3] shows results for various values of /is; qb and 
bs are the same sign for the values on the left, and the 
opposite sign for the values on the right. The profile that 
is closest to being current-free is /ts =0, and indeed we 



find there that even for Ha=200 there is no sign of any 
destabilizing infiuence of the field, for either axisymmet- 
ric or nonaxisymmetric perturbations. For < fls < Ai 
the magnetic field stabilizes the flow for both m = Q and 
m — 1. 

For jli < fiB < Ao the m — 1 mode should be unstable, 
while the m = mode should be stable. The values 
fiB — ^ and [jLB = are examples of this situation. There 
is always a crossover point at which the most unstable 
mode changes from m = to m = 1. Note also how for 
/is = 1, the critical Reynolds number increases for the 
TO = mode, before suddenly decreasing for the to = 1 
mode (Fig. [3l left bottom plot). We have the interesting 
situation therefore that weak fields initially stabilize the 
TC-fiow, before stronger fields eventually destabilize it, 
via a non-axisymmetric mode. Beyond Ha = 150 (the 
same value we saw previously in Fig. [2]), the flow is 
unstable even for Re — 0. 

Except for the almost current-free profile fiB — 0, all 
other values share this feature, that there is a critical 
Hartmann number beyond which the basic state is un- 
stable even for Re = 0. Let Ha^"^ and Ha*^^-* denote these 
critical Hartmann numbers, for m = and 1 respectively. 
For the profiles with the largest gradients both modes are 
unstable. Strikingly, in these cases m = is always more 
unstable than to = 1, that is, Ha^"^ < Ha^^^ see the plots 
for fiB = 4: and fiB = —2 of Fig. [31 



B. The required electric currents 

Let /axis be the axial current inside the inner cylinder 
and /fluid the axial current through the fluid (i.e. between 
inner and outer cylinder). The toroidal field amplitudes 
at the inner and outer cylinders are then 



-Bin — 



5i?in' 



J-, (/axis + /fluid) /'no^ 

Bout - =^ , (,2dj 

Oitout 



where /?, B and / are measured in cm. Gauss and Am- 
pere. Expressing /axis and /fluid in terms of our dimen- 
sionless parameters one finds 



5Ha 



1/2 



1/2 



and 



/fluid — 



(24) 



(25) 



in Ampere. Note also how the required currents depend 
on the radius ratio f), but not on the actual physical di- 
mensions Rin and Rout- Making the entire device bigger 
thus reduces the current density, inversely proportional 
to the square of the size. By making the device suf- 
ficiently large one can thereby prevent Ohmic heating 
within the fluid from becoming excessive. 

The results for the critical Hartmann numbers are now 
applied to two different conducting liquid metals, sodium 
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FIG. 2: Marginal stability curves for m = (dashed) and m — 1 (solid). The hatched area is thus the region that is stable to 
both. Magnetic Prandtl numbers as indicated on each plot. 



TABLE I: Material parameters of liquid metals that might be 
used for magnetic TC-experiments. 





P 


[g/cm^] 


u [cm^/s] 


Tj [cm^/s] 




sodium 




0.92 


7.10 ■ 10"=* 


0.81 ■ 10^ 


8.15 


gallium-indium-tin 




6.36 


3.40 ■ 10"^ 


2.43 ■ 10^ 


25.7 



and gallium- indium-tin [3], whose material parameters 
are given in Table [J We also wish to consider the effect 
of varying the radius ratio i). Tables Ultil VI give the values 
of the electric currents needed to reach the lesser of Ha*^"^ 
and Ha(^\ for the three values i) = 0.25, 0.5 and 0.75, and 
for jlB ranging from —10 to 10 in each case. Note that 
for large l/fsl, Ha^*^^ scales as l/As, and /fluid approaches 
a constant value. The calculated currents are lower for 
fluids with smaller y^JlQpUrj (i.e. sodium is better than 
gallium) . 

In these Tables, the most interesting experiment, with 
the almost uniform field fiB = 1 (see Fig. [31 top-right) 
is indicated in bold. For a container with a medium gap 
oi fj = 0.5, parallel currents along the axis and through 
the fluid of 6.16 kA for sodium and 19.4 kA for gallium 
are necessary. Such sodium experiments should indeed 
be possible. Experiments with a wider f} — 0.25 gap are 
even easier; in that case even gallium experiments should 
be possible, with a current of 9.29 kA required (see Table 

mi). 



The Reynolds numbers that would be required to ob- 
tain not just these Re — pure Tayler instabilities, but 
also the transition points from m = to to = 1 are also 
not difficult to achieve; for i?out ^ 10 cm, say, rotation 
rates of order 10~^ Hz are already enough. 



TABLE II: Characteristic Hartmann numbers and electric 
currents for a wide gap container (?) = 0.25) with conducting 
walls, using either sodium or gallium-indium-tin (in brackets). 



As 


Ha(0) 


Ha(^) 


/axis [kA] 


/fluid 


[kA] 


-10 


2.29 


2.05 


0.0483 (0.152) 


-1.98 


(-6.24) 


-5 


4.23 


3.98 


0.0937 (0.296) 


-1.97 


(-6.21) 


-4 


5.14 


4.93 


0.116 (0.366) 


-1.97 


(-6.22) 


-3 


6.63 


6.50 


0.153 (0.483) 


-1.99 


(-6.27) 


-2 


9.69 


9.71 


0.228 (0.721) 


-2.05 


(-6.49) 


-1 


24.7 


22.5 


0.530 (1.67) 


-2.65 


(-8.35) 


1 


oo 


41.7 


0.982 (3.10) 


2.94 


(9.29) 


2 


oo 


13.8 


0.325 (1.02) 


2.27 


(7.17) 


3 


oo 


8.27 


0.195 (0.614) 


2.14 


(6.75) 


4 


oo 


5.93 


0.140 (0.440) 


2.09 


(6.60) 


5 


10.8 


4.63 


0.109 (0.344) 


2.07 


(6.53) 


10 


3.23 


2.205 


0.0519 (0.164) 


2.02 


(6.39) 
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FIG. 3: The marginal stability curves for m — (dashed) and m — 1 (solid). Pm= 10~^, f) = 0.5, /in ~ 0, and /is as indicated. 
Note also how the critical Reynolds numbers are always of the same order of magnitude as the nonmagnetic result 68, which 
is very easy to achieve in the laboratory. 



TABLE III: The same as in Table|lT]but for r) = 0.5 TABLE IV: The same as in Table|lT]but for i) = 0.75 (narrow 

gap container). 



/is 


Ha(o' 


Ha(^) 


/axis [kA] 


/fluid [kA] 


-10 


3.96 


5.02 


0.161 (0.509) 


-3.39 (-10.7) 


-5 


7.73 


9.85 


0.315 (0.994) 


-3.47 (-10.9) 


-4 


9.61 


12 


0.392 (1.24) 


-3.53 (-11.1) 


-3 


12.8 


16.2 


0.522 (1.65) 


-3.65 (-11.5) 


-2 


19.8 


24.8 


0.807 (2.55) 


-4.04 (-12.7) 


-1 


59.3 


63.7 


2.42 (7.63) 


-7.25 (-22.9) 


1 


oo 


151 


6.16 (19.4) 


6.16 (19.4) 


2 


oo 


35.3 


1.44 (4.54) 


4.32 (13.6) 


3 


21.0 


20.6 


0.840 (2.65) 


4.20 (13.2) 


4 


13.2 


14.6 


0.538 (1.70) 


3.77 (11.9) 


5 


9.84 


11.4 


0.401 (1.27) 


3.61 (11.4) 


10 


4.44 


5.4 


0.181 (0.571) 


3.44 (10.8) 



As 


Ha(«) 


Ha(i) 


/axis [kA] 


/fluid [kA] 


-10 


9.27 


12.2 


0.655 (2.06) 


-9.39 (-29.6) 


-5 


18.6 


24.1 


1.31 (4.13) 


-10.0 (-31.5) 


-4 


23.3 


30.1 


1.65 (5.20) 


-10.5 (-33.1) 


-3 


31.6 


40.3 


2.23 (7.03) 


-11.1 (-35.0) 


-2 


50.4 


63.4 


3.56 (11.2) 


-13.1 (-41.3) 


-1 


163. 


177. 


11.5 (36.3) 


-26.8 (-84.5) 


1 


oo 


632 


44.6 (141) 


14.9 (47.0) 


2 


66.8 


87.3 


4.72 (14.9) 


7.87 (24.8) 


3 


36.5 


49.8 


2.58 (8.14) 


7.74 (24.4) 


4 


25.8 


35.2 


1.82 (7.89) 


5.74 (24.9) 


5 


20.1 


27.3 


1.42 (4.48) 


8.05 (25.4) 


10 


9.63 


13.0 


0.680 (2.14) 


8.39 (26.5) 



C. Insulating cylinder walls 



Calculations were also done for insulating cylinder 
walls; the results are given in Fig. [D They are gen- 
erally similar as those for the conducting cylinders, but 
with one important exception. The m = and m = 1 
stability curves now almost always cross each other, as 
they do for conducting cylinders only for almost current- 



free i?0 profiles. One can again observe for not too steep 
profiles (for fiB — ±1) how for weak fields the m — 
mode stabilizes the rotation until beyond the cross-over 
point the m = 1 mode strongly destabilizes the rotation. 
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FIG. 4: The same as in Fig. |3]but for insulating cylinder walls. 



TABLE V: The same as in Table |TT] (wide gap, rj = 0.25) but TABLE VL The same as in Table |TTT] (medium gap, ?) = 0.5) 
for insulating cylinders. but for insulating cylinders. 



flB 


Hao 


Hai 


/axis [kA] 


/fluid 


[kA] 




Hao 


Hai 


/axis [kA] 


/fluid 


[kA] 


-10 


3.63 


1.79 


0.0421 (0.133) 


-1.73 


(-5.45) 


-10 


6.09 


4.66 


0.190 (0.599) 


-3.99 


[-12.6) 


-5 


6.65 


3.55 


0.0836 (0.264) 


-1.76 


(-5.54) 


-5 


11.8 


9.31 


0.380( 1.20) 


-4.18 


[-13.2) 


-4 


8.04 


4.42 


0.104 (0.328) 


-1.77 


(-5.58) 


-4 


14.6 


11.7 


0.477 (1.50) 


-4.29 


[-13.5) 


-3 


10.3 


5.89 


0.139 (0.437) 


-1.81 


(-5.68) 


-3 


19.4 


15.7 


0.640 (2.02) 


-4.48 


[-14.1) 


-2 


14.7 


8.91 


0.210 (0.662) 


-1.89 


(-5.60) 


-2 


29.3 


25.2 


1.03 (3.24) 


-5.15 


[-16.2) 


-1 


30.6 


20.5 


0.483 (1.52) 


-2.42 


(-7.60) 


-1 


73.6 


64.6 


2.63 (8.31) 


-7.89 


[-24.9) 


1 


oo 


30.7 


0.723 (2.28) 


2.17 


(6.84) 


1 


oo 


109. 


4.44 (14.0) 


4.44 


[14.0) 


2 


oo 


10.7 


0.252 (0.794) 


1.76 


(5.56) 


2 


oo 


28.1 


1.15 (3.61) 


3.45 


[10.8) 


3 


oo 


6.63 


0.156 (0.492) 


1.72 


(5.41) 


3 


32.6 


17.2 


0.701 (2.21) 


3.51 


[11.1) 


4 


oo 


4.83 


0.114 (0.359) 


1.71 


(5.39) 


4 


20.5 


12.5 


0.510 (1.61) 


3.57 


[11.3) 


5 


17.4 


3.81 


0.0897 (0.283) 


1.70 


(5.38) 


5 


15.3 


9.81 


0.400 (1.26) 


3.60 


[11.3) 


10 


5.18 


1.86 


0.0438 (0.138) 


1.71 


(5.38) 


10 


6.86 


4.78 


0.195 (0.615) 


3.71 


[11.7) 



V. CONCLUSIONS 

We have shown how complex the interaction of mag- 
netic fields and differential rotation can be in Taylor- 
Couette flows, including also a strong dependence on the 
magnetic Prandtl number. For large Pm the field desta- 
bilizes the differential rotation, whereas for small Pm it 
stabilizes it. However, if the field (or rather the cur- 
rent) is too great, then the Tayler instabilities will always 



destabilize any differential rotation. 

In order to prepare laboratory experiments, we also did 
calculations at values of Pm appropriate for liquid met- 
als, for both conducting and insulating cylinder walls. In 
particular, we considered the almost uniform field profile 
fiB — ^- For both conducting and insulating boundaries, 
the field is initially stabilizing, but after the most unsta- 
ble mode switches from m = to m = 1 it is strongly 
destabilizing, until the pure Tayler instability sets in even 
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TABLE VII: The same as in Table IWl (narrow gap, i) = 0.75) 
but for insulating cylinders. 





Mao 


riai 


J axis [KAJ 


r A 1 
J fluid [kaj 


-lU 


1/10 


lo.o 


n 07r^ Q (Y7\ 

U.y / (^o.U / j 


-14. U (^-44. UJ 


c 
-0 


02 


07 7 


1 1 '-7\ 


1 n /1 7 


A 

-4 


00.4 


'XA Q 


/Ifi /'7 77^ 
Z.4D \i ' i ( ) 


1 /' AO 0^ 


o 
-o 




/17 
4 ( .Z 


O.oo ^^1U.0J 


-iO. / l^-Oz.O J 


o 
-z 


1A ft 


7/1 7 
( 4. / 






-1 


208. 


205. 


14.5 (45.7) 


-33.8 (-107) 


1 


oo 


464. 


32.8 (103.) 


10.9 (34.3) 


2 


103. 


80.3 


5.67 (17.9) 


9.45 (29.8) 


3 


56.2 


49.1 


3.47 (10.9) 


10.4 (32.7) 


4 


39.7 


35.9 


2.54 (7.99) 


11.0 (34.6) 


5 


30.9 


28.4 


2.18 (6.32) 


12.4 (35.8) 


10 


14.8 


13.9 


0.982 (3.10) 


12.1 (38.2) 



at Re = 0. 



For various gap widths and field profiles, we also com- 
puted the critical Hartmann numbers and the corre- 
sponding electric currents. Tables II-VII give the re- 
quired currents for both conducting and insulating walls; 
note how insulating walls (Tables V-VII) typically re- 
quire lower currents than conducting walls (Tables II- 
IV). The other clear trend is that the currents are lesser 
for wider gaps and greater for narrower gaps. An opti- 
mal experiment might therefore have rj = 0.25, insulating 
walls, and fiB = 1, which would require only 6.84 kA even 
with gallium- indium-tin (Table V). 



Acknowledgments 

The work of DS was supported by the Helmholtz In- 
stitute for Supercomputational Physics in Potsdam and 
partly by the Leibniz Gemeinschaft under program SAW. 



[1] J. A. Sellwood and S. A. Balbus, Astrophys. J. 511, 660 
(1999). 

[2] N. Dziourkevitch, D. Elstner, and G. Riidiger, Astron. 

Astrophys. 423, L29 (2004). 
[3] F. Stefani, T. Gundrum, G. Gerbeth, G. Riidiger, M. 

Schultz, J. Szklarski, and R. HoUerbach, Phys. Rev. Lett. 

97, 184502 (2006). 
[4] G. Riidiger, R. HoUerbach, F. Stefani, T. Gundrum, G. 

Gerbeth, and R. Rosner, Astrophys. J. 649, L145 (2006). 
[5] R. HoUerbach and G. Riidiger, Phys. Rev. Lett. 95, 

124501 (2005). 

[6] R. HoUerbach, G. Riidiger, M. Schultz and D. Elstner, 

Phys. Rev. Lett, (submitted). 
[7] R. J. Tayler, Mon. Not. R. Astron. Soc. 161, 365 (1973). 
[8] H. C. Spruit, Astron. Astrophys. 349, 189 (1999). 
[9] N. A. Krivova and S. K. Solanki, Astron. Astrophys. 394, 



701 (2002). 

[10] H. Korhonen, S. V. Berdyugina, K. G. Strassmeier, and 
I. Tuominen, Astron. Astrophys. 379, L30 (2001). 

[11] J. Braithwaite and A. Nordlund, Astron. Astrophys. 450, 
1077 (2006). 

[12] J. Braithwaite, Astron. Astrophys. 453, 687 (2006). 

[13] D. Michael, Mathematika 1, 45, (1954). 

[14] S. Chandrasekhar, Hydrodynamic and Hydromagnetic 

Stability (Clarendon, Oxford 1961). 
[15] E. P. Velikhov, Sov. Phys. JETP 9, 995 (1959). 
[16] G. Riidiger and D. Shalybkov, Phys. Rev. E 66, 016307 

(2002). 

[17] P. A. Gilman and P. A. Fox, Astrophys. J. 484, 439 
(1997). 



